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find a consistent explanation for both recently measured CP asymmetry in 



Qh. B° — * 4>Ks decay and the large branching ratio of — * rj' decay, which 



are inconsistent with the Standard Model (SM) prediction. We also investi- 



gate other charmless hadronic B — > PP and B — > VP decay modes whose 

u : 

experimental data favor the SM: for instance, recently measured CP asym- 
metries in B° — ► n'Ks and £? — > J/^/Kg. We find that all the observed data 
can be accommodated for certain values of Jjip couplings. 
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The main mission of B factories, such as KEK— B and SLAC— B, is to test the Standard 
Model (SM) and further to find possible new physics effects in B meson systems. For 
this goal, a variety of useful observables can be measured to be compared with theoretical 
expectations. One of the important observables is CP asymmetries in various B meson 
decays, and another is branching ratios (BRs) for rare B decay processes. 

CP violation in B system has been confirmed in measurements of time-dependent CP 
asymmetries in B — > J/^>K S decay. The world average of the asymmetry in B — > J/^>K S 
[1] is given by 

sin(2/3) j/y Ks = 0.734 ± 0.054 , (1) 

which is consistent with the SM expectation. However, recent measurements of sin(2/3) in 
B° — > (f)Ks decay disagree with the above value. The measurements of time-dependent 
CP asymmetries in B° — > <pK s have been recently reported by BaBar [2] and Belle [3], 
respectively: 

sin(2/3)gg- = -0.19l8;ijj ±0.09, (2) 
sin(2/3) = -0.73 ± 0.64 ± 0.18 . (3) 

Both experimental results show (a tendency of) negative values of sin(2/3). 

Because within the SM the difference between the asymmetries in B — > J/^>K S and 
B° — > 4>K S is expected to be [4] 

| sin(2/3) j/^ Ks - sin(2/5)^ s | < 0(A 2 ) (4) 

with A 0.2, these results indicate 2.7 a deviation from the SM prediction and may reveal 
new physics effects. 

The decay process B° — > is a pure penguin process b — > sss, and in the SM the 
time-dependent CP asymmetry of this mode is expected to measure the same sin(2/3) as in 
B — > J/^/Ks, as shown in Eq. (4). If there are any new physics effects in the measurements 
of time-dependent CP asymmetries, the effects can arise from new contributions to the 
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B° — B° mixing amplitude and/or the decay amplitude of each mode. Because the new 
physics effects to the B° — B° mixing would be universal to all the neutral B® meson decays, 
the discrepancy between the measured values of sin(2/3) j/^k s and sin(2/5)^ s may indicate 
the possibility of new physics effects in the decay amplitude of B° — > 4>K S . Supposing that 
the mode B° — > 4>K S reveals new physics effects [4,5], one might expect that other modes 
having the same internal quark level process b — > sss (e.g., B° — > r/'Ks) would reveal the 
similar effects. Interestingly, the recent measurements of CP asymmetries in B° — > rj'Ks 
by Belle [3] agree well with the results of sin(2/9) j/^>k s given in Eq. (I) (see Table I). 
Therefore, any successful explanation invoking new physics should accommodate the data 
disfavoring the SM, such as sin(2/9)^ g , and simultaneously all the data favoring the SM, 
such as sin(2 / 9)^i S : s and sin(2/9) j/^k s - Motivated by the recent measurements of sin(2/5)^ s , 
a few first analyses [1,6-9] have been done. But, all of them explain only the deviation of 
sin(2/5)0x s from the SM prediction, or in addition, make (at most) qualitative statements 
on sm(2f3) v < Ks [1,6-9]. 

There also exist plenty of experimental data observed for charmless hadronic B — > PP 
(P denotes a pseudoscalar meson) and B — > VP (V denotes a vector meson) decay modes, 
which are well understood within the SM. However, among the B — > PP decay modes, the 
BR of the decay mode B^ — > 7]'K ± is found to be still larger than that expected within the 
SM 1 [10-13]: the experimental world average is B{B ± -> r/K*) = (75 ± 7) x 1(T 6 . The 
SM contribution is about 3<x smaller than the experimental world average [13]. Among the 
B — > VP decay modes, the experimentally observed BR for the decay B° — > 7]K*° has been 
also found to be 2a larger than the SM expectation [13]. 

In this letter we try to find a consistent explanation for all the observed data in charmless 
hadronic B — > PP and B — > VP decays in the framework of -R-parity violating ($ p ) 

1 Attempts to explain the large B{B ± — >■ n'K^) within the SM have been made: e.g., by using the 
anomalous g — g — n' coupling [14]. 
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super symmetry (SUSY). In ft p SUSY, the effects of ft p couplings on B decays can appear at 
tree level and can be in some cases comparable to the SM contribution. Our main focus is on 
the recent measurement of CP asymmetry in B° — > <pKs and the large BR for B ± — > r)'K ± : 
both data appear to be inconsistent with the SM prediction. In order to achieve this goal, 
we investigate all the observed B — > PP and B — > VP decay modes and explicitly show that 
indeed all the observed data can be accommodated in a consistent way for certain values of 
ftp couplings. 

The ftp part of the superpotential of the minimal supersymmetric standard model 
(MSSM) can contain terms of the form 

}% = n t L t H 2 + \ l]k UL 3 El + K jk L t Q,D c k + X^U'Dpi (5) 

where E iy Ui and D { are respectively the i-th type of lepton, up-quark and down-quark 
singlet superfields, Lj and Qi are the SU(2)l doublet lepton and quark superfields, and H 2 
is the Higgs doublet with the appropriate hypercharge. From the symmetry reason, we need 
\jk = —Xjik and \"j k = —\" k y The bilinear terms can be rotated away with redefinition of 
lepton and Higgs superfields, but the effect reappears as As, A's and lepton-number violating 
soft terms [15]. The first three terms of Eq. (5) violate the lepton number, whereas the 
fourth term violates the baryon number. We do not want all these terms to be present 
simultaneously due to catastrophic rates for proton decay. In order to prevent proton decay, 
one set needs to be forbidden. 

For our purpose, we will assume only A'— type couplings to be present. Then, the effective 
Hamiltonian for charmless hadronic B decay can be written as [16] 

Heff(b ~ y djdkdn) = df kn [d na ^d j/3 d kj3 ^^ R b a \ + df kn [d na ^b/3 d k ^^ R d ja ] , 

H eff( b ~^ UjU k d n ) = uf kn [u ka ^U j/3 dnp'J^Rba} , (6) 

with 
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3 y y* 3 y y* 

A ijk A in3 iL \" A i3k A inj 



^ = EW «.*.» = 1.2) 



I/, 



3 y y* 

_ \ ^ A ijn ik3 



ElS^ (** = l,n = 2), 



i=l C "' i e iL 

where a and /? are color indices and 7^ L = 7 M (1 ± 75). The leading order QCD correction 
to this operator is given by a scaling factor / ~ 2 for = 200 GeV. 

The available data on low energy processes can be used to impose rather strict constraints 
on many of these couplings [17-19]. Most such bounds have been calculated under the 
assumption of there being only one non-zero ft p coupling. There is no strong argument to 
have only one ftp coupling being nonzero. In fact, a hierarchy of couplings may be naturally 
obtained [18] on account of the mixings in either of the quark and squark sectors. In this 
work, we try to find out the values of such ft p couplings for which all available data are 
simultaneously satisfied. An important role will be played by the \' 32i -type couplings, the 
constraints on which are relatively weak. The BR of B — > X s uu can put bound on A^A^ 
in certain limits. Using Ref. [20] and the experimental limit (BR< 6.4 x 10~ 4 ) on the BR 
of B — > X s vv [21], we find that A' < 0.07 (for = 200 GeV). However, if we go to any 
realistic scenario, for example grand unified models (with _R-parity violation), we find a 
natural hierarchy among the sneutrino and squark masses. The squark masses are much 
heavier than the sneutrino masses and the bound does not apply any more for rxiy = 200 
GeV. 

In our calculation, we use the effective Hamiltonian and the effective Wilson coefficients 
(WCs) for the processes b — > sqq' and b — > dqq' given in Ref. [22]. The decay amplitude of 
B~ — > <fiK~ decay mode is given by 

A^k = A^k + Af v K , (8) 

where the SM part of the amplitude is 

G 111 

= -^=VtbV t * s (a 3 + a 4 + a 5 - -a 7 - -a 9 - -a 10 )A^ , (9) 

and the ftp part of the amplitude involves only 0^22 an d ^22 : 



A% = (d$ 22 + dg 2 ) [ZA+] , (10) 

with A<p = (K\s^(l — 7 5 )6|S) (0|s7 M s|O). This particular structure of A<p is obtained from 
the operators d no rf1m)djp dkpl nR(L)b a given in Eq. (6), which are derived from the operators 
s l(r)S r(l)S R(L)bL(R) by Fierz transformation [23]. Here the effective coefficients Oj are defined 
as a,i = c e /* + £cf+{ (i — 0( id) and a, = c e /* + £cf^{ (i = even) with the effective WCs 
at the scale m b [22], and by treating £ = l/iV c (iV c denotes the effective number of color) as 
an adjustable parameter. 

The CP asymmetry parameter sm(2/3) XY for B — > XY can be written as 

S m(2/3W = - n 2 '"' A - Y " , (11) 
(1 + \\xy\ z ) 

where A xy = e~ 2i ^ {A XY / A XY ) = e~ i( - 2f3+ ^ \A XY /A XY I with i iy /i xy = e - i9 |^ xy /^l X y|. 
Here (3 denotes the CP angle in the SM and (3 denotes the effective CP angle given by 
2/3 = 2/3 + 9 with a new weak phase 6 arising from new physics. A X y = A x ^ + A X P Y and 
A X y is its CP-conjugate amplitude. For B° — > (bKs, in the SM, A$k s I A^Ks — 1, so (3 — (3. 
But, in % SUSY, generally |A?>k s /A^ s | ^ 1 and 6» 7^ 0. 
The fc p part of the amplitude of B~ — > decay is 

= (^121 ~~ ^112) + (^222 ~~ ^222) ~ (-^v ~~ ~~ i^ri' 



4. fr/ L _ r/ R ^ —4" + 7, R 



(m s + m u )(m b - m u ) 



(12) 



where m = m 2 ,/ / (m& — m s ) and 



A u(s) = ff) FB ^ K{m 2 B _ rn 2 ch 

A K = f K F B ^(m 2 B -rr%). 
pB^K(ri') (-[gno^gg the hadronic form factor for B — > K(rj') and /k"(t/') is the decay constant 
of K(r)') meson 2 . Analogous expressions hold for B ± — > r\K ± where we have to replace A™, 



2 The definition of /^/^ is given in Ref. [22], considering the r] — rj mixing. The mixing angle 
9 = -22° is used. 
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TABLE I. CP asymmetries in the decay modes B° — > 4>Kg and B° — > r/ifs. 



sin(2/3) 


Case 1 


Case 2 


experimental data 


sin(2/3)^ s 
sin(2^),,/^ s 



0.73 


-0.82 
0.72 


— 0.19±g;|g ± 0.09 (BaBar), -0.73 ± 0.64 ± 0.18 (Belle) 
0.76 ± 0.36±g;g| (Belle) 



by A^, by and m v > by m v . Replacing a pseudoscalar meson by a vector meson, we 
also get similar expressions for the amplitudes of — > rj'K*^ ^ modes. 

In Eqs. (10) and (12), we note that A^ K includes the (d 222 + ^222) term only, while 
^ includes the {d L 222 - dg 2 ) term. Furthermore, in the % SUSY framework, we find that 
the positive values of d 222 and negative values of d\ 22 can increase the BR for the process 
_g± V [K ± , keeping most of the other B — > PP and 5 — > VP modes unaffected. The 
other ^ p combinations are either not enough to increase in the BR for B^ — > r]'K ± or affect 
too many other modes. From now on, we will concentrate on contributions of d 222 an d ^222 
which are relevant to the process b — > sss. We divide our results into the following two 
cases. 

Case 1: We use the following input parameters: the CP angle 7 = 110°, the s quark 
mass m s (at m h scale) = 85 MeV, and the decay constants and the hadronic form factors 
used in Ref. [13]. We take (5 = 26° as our input, which corresponds to sin(2/?)j/$x = 0.78 
given in Eq. (1). We set <i 222 = ke~ t9 ' and d 222 = —ke t9> , where k = \d 222 \ = \d 222 \ 
and 9' is a new weak phase arising from d 222 oc A- 32 A- 22 , or (d 222 )* oc A- 22 A- 23 . In this 
case, A$x is purely imaginary and introduces a new weak phase to the decay amplitude of 
B — > (j)K modes, which can cause non-zero direct CP asymmetries, different from the SM 
prediction. In contrast, A^ K introduces no new phase in this choice of d 222 and d 222 . We 
find a constructive contribution to the SM part of the amplitude for B + — > i]'K + which 
helps to satisfy the experimental data. 

For |A 322 | = |A 332 | = |A 323 | = 0.055, tan 9' = 0.52, and m susy = 200 GeV, we find 
sin(2/3)^ g = and sin (2/?)^^ = 0.73 for £ = 0.45, as shown in Table I. This result agrees 
well with the recently reported experimental data. Figures 1 and 2 show the BR versus £ for 
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B + — > 0K + and 5 + — > i]'K + , respectively. In both figures, the shaded region represents the 
allowed region by the experimental results and the solid line corresponds the SM prediction. 
The SM BR of B + — > 7]'K + is far below the experimental lower bound. Our result in 
this case is denoted by the dotted line, which is clearly consistent with the experimental 
data: B{B + -> r]'K + ) = 69.3 x 1CT 6 for f = 0.45. In Table II, we estimate the BRs and 
CP (rate) asymmetries Acp for B — > rf'^K^ and B + — > (f)K + , where Acp is defined by 
Acp = w hh b and / denoting b quark and a generic final state, respectively. 

The estimated BRs are well within the experimental limits. The BR of B° — > (pK° has 
also been measured, B{B° -> <f)K°) = (8.1^ ± 0.8) x 10~ 6 [2]. The theoretical value for 
the BR of this mode is same as that of B + — > (j)K + and is allowed by the experimental 
data. Further, We note that Acp's for B + — > <j)K + and 5° — > ^fT* are large and have 
positive signs. This can be taken as a prediction of this model. It is hopeful that the CDF 
at Fermilab could measure the CP asymmetry of B + — > (j)K + mode using their two track 
trigger in the ongoing RUN II experiment [24]. 

So far we have assumed that the magnitudes and phases are same for d 222 and d 222 . 
We can assume them to be different and still obtain good fits. For example, if we set 
^222 — k L e~ ie L and df 22 = -k R e id> R and choose: k L = 1.0 x 10~ 8 , k R = 8.0 x 10~ 9 and 
tan^ = 0.4, tan 0^ = 0.52, we find sin(2/^) 0X5 = 0.017 and B(B ± -> QK*) = 7 x 10- 6 . 
The BRs of B — > rj'K modes remain unchanged. 

The other observed B — > PP and B — > VP decay modes without or in the final 
state, such as B — > 7r7r, 5 — > — > p7r, 5 — > c<j7r, and so on, are not affected in this 
scenario, because 0^22 an d ^222 are relevant to the process b — > sss only. The estimated BRs 
for those modes by using the above input values are consistent with the experimental data 
for f = 0.45 [13]. 

Case 2: In Case 1, we have generated the very small value for sin(2/5)0x s - Let us now try 
to generate a large negative sin(2/9)0x s , since one experiment is preferring a small negative 
value and the other experiment is preferring large. We use the smaller value of 7 and m s : 
7 = 80° and m s (at m& scale) = 75 MeV, keeping the other inputs unchanged. We also use 
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FIG. 1. BR vs £. The dotted and dot-dashed lines correspond to Case 1 and Case 2 respectively. 
The solid line corresponds to the SM (The SM BR is same for both cases). The shaded region is 
allowed by the experimental data. 



TABLE II. The branching ratios (B) and CP rate asymmetries (Acp) for B — > rj^'K^ and 
B -> 





Case 1 




Case 2 




mode 


B x 10 6 




B x 10 6 


A C p 




69.3 


-0.01 


76.1 


-0.01 


-> ^ *+ 


27.9 


-0.04 


35.2 


-0.03 




107.4 


0.00 


98.9 


0.00 


5° -> 7]K*° 


20.5 


0.71 


11.7 


0.15 


B + 


8.99 


-0.21 


8.52 


-0.25 
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7.5 



b 

i — i 

+ £ 

T 

+i 

PQ 



« 2.5 








FIG. 2. BR vs £. The dotted and dot-dashed lines correspond to Case 1 and Case 2 respectively. 
The solid lines correspond to the SM. The upper solid line is for Case 1 and the lower solid line is 
for Case 2. The shaded region is allowed by the experimental data. 
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the same form of d\ 22 and as in Case 1: the 1/lp part of the decay amplitude introduces 
a new weak phase in the mode B — > <pK, but no phase in the process B + — > i]'K + . 

In this case, we find that a consistent explanation of all the observed data for B — > PP 
and 5 — > VP is possible for |A3 22 | = (A^l = | Ag 23 1 = 0.069, tan#' = 2.8, and m susy = 
200 GeV. In this scenario, a large value of sin (2/3) with the negative sign is possible for 
B° -> <j>K 3 \ sin(2^)^ g = -0.82 for f = 0.25, as shown in Table I. The value of sin(2/3) for 
B° -> is similar to that of Case 1: sin(2^)^ Ks = 0.72 for f = 0.25. The estimated 

BRs for B + — > and 5 + — > ?/fr + are denoted by the dot-dashed lines in Figs. 1 and 2, 
respectively: B(B + -> 0fsT + ) = 8.52 x 10" 6 and £(£+ -> ?/^ + ) = 76.1 x 10~ 6 for £ = 0.25. 
The BRs of other _B — > PP, VP modes satisfy the experimental data for £ = 0.25 [13]. 

In conclusion, we have shown that in R p violating SUSY it is possible to understand 
consistently both recently measured CP asymmetry sin(2/3) in B° — > <pK s decay and the 
large branching ratio of B ± — > r/'K* 1 decay, which appear to be inconsistent with the SM 
prediction. Experimental results imply that CP asymmetry for B — > <pK s is different from 
the SM value, but for B — > i]'K s it is same as in the SM. We have found that this tendency 
can be consistently understood in the framework of this model. We have searched for possible 
parameter space and found that all the observed data for B — > PP and B — > VP decays 
can be accommodated for certain values of Ift p couplings. Experimentally, these couplings 
can be observed through the direct production of SUSY particles at Tevatron (e.g., xfxl 
production and their subsequent decays to final states with multiple b quarks) and other 
rare decays like B s — > /i ± / u T [25]. 
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